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ZEROTREE ENCODING OF WAVELET DATA 
BACKGROUND 

The invention generally relates to zerotree encoding of wavelet data, such as 
zerotree encoding of wavelet coefficients, for example. 

Data compression typically removes redundant information from a set of data to 
produce another set of data having a smaller size. This smaller size may be beneficial, 
for example, for purposes of transmitting the data over a bus or network. 

For example, the pixel intensities of an image may be indicated by a set of 
coefficients, and these coefficients may be represented by digital image data. For 
purposes of compressing the image data, the data may be transformed to reveal redundant 
information, i.e., information may be removed via data compression. For example, the 
image data may be transformed pursuant to a wavelet transformation, a transformation 
that effectively decomposes the image into spatially filtered images called frequency 
subbands. In this manner, the subbands may reveal a significant amount of redundant 
information that may be removed by compression techniques. 

Referring to Fig. 1, as an example, image data that indicates pixel intensities of an 
image 12 may undergo wavelet transformations to decompose the image 12 into 
subbands. Due to the nature of the transformations, the subbands appear in different 
decomposition levels (levels 14, 16 and 18, as examples). In this manner, to decompose 
the original image 12 into subbands 14a, 14b, 14c and 14d of the first decomposition 
level 14, the one dimensional Discrete Wavelet Transform (DWT) is applied row-wise 
and then column-wise. In one dimensional DWT, the signal (say a row-wise) is first low- 
pass filtered and sub-sampled by dropping the alternate filtered output to produce the . 
low-frequency subband (L) which is half the size of the original signal. Then the same 
signal is high-pass filtered and similarly sub-sampled to produce the high-frequency 
subband (H) which is half the size of the original signal. When the same one dimensional 
operation is applied column-wise on the L subband, it produces two subbands LL and 
LH. Similarly, applying the same one dimensional operation column-wise on the H 
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subband, it produces two subbands HL and HH subbands. As a result after two- 
dimensional Discrete Wavelet Transform, the original image 12 is decomposed into four 
subbands: the LL subband 14a 5 the LH subband 14b, HL subband 14c and HH subband 
14d. Sizes of the row and column of each of these subbands is half the sizes of the row 
5 and column of the original images due to the sub-sampling operation. The values of these 
subbands are called the wavelet coefficients and hence the subbands may be represented 
by an associated matrix of wavelet coefficients. 

The LL subband 14a indicates low frequency information in both the horizontal 
and vertical directions of the image 12 and typically represents a considerable amount of 
10 information present in the image 12 because it is nothing but the sub-sampled version of 
the original image 12. The LH subband 14b indicates low frequency information in the 
horizontal direction and high frequency information in the vertical direction, i.e., 
% i horizontal edge information. The HL subband 14c indicates high frequency information 

W in the horizontal direction and low frequency information in the vertical direction, i.e., 

□ 15 vertical edge information. The HH subband 14b indicates high frequency information in 

|S the horizontal direction and high frequency information in the vertical direction, i.e., 

u 1 diagonal edge information. 

P Since LL subband 14a is nothing but the sub-sampled version of the original 

tt image, it maintains the spatial characteristics of the original image. As a result, the same 

20 DWT decomposition can be further applied to produce four subbands that have half the 
resolution of the LL subband 14a in both the vertical and horizontal directions: the LL 
subband 16a, LH subband 16b, HL subband 16c and HH subband 16d. Hence the LL 
subband 16a is again the sub-sampled version of the LL subband 14a. Hence LL subband 
16a can be further decomposed to four subbands that have half of its resolution in both 
25 horizontal and vertical directions: LL subband 18a, LH subband 18b, HL subband 18c 
and HH subband 18d. 

The subbands of the lower decomposition levels indicate the information that is 
present in the original image 12 in finer detail (i.e., the subbands indicate a higher 
resolution version of the image 12) than the corresponding subbands of the higher 
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decomposition levels. For example, the HH subband 18d (the parent of the HH subband 
16d) indicates the information that is present in the original image 12 in coarser detail 
than the HH subband 16d (the child of the HH subband 18d), and the HH subband image 
14d (another descendant of the HH subband 18d) indicates the information that is present 
5 in the original image 12 in finer detail than the HH 16d and 18d subbands. In this 
manner, a pixel location 24 of the HH subband image 18d corresponds to four pixel 
locations 22 of the HH subband 16d and sixteen pixel locations 20 of the HH subband 
14d. 

Due to the relationship of the pixel locations between the parent subband and its 

10 descendants, a technique called zerotree coding may be used to identify wavelet 

coefficients called zerotree roots. In general, a zerotree root is a wavelet coefficient that 
satisfies two properties: the coefficient has an insignificant intensity, and all of the 
descendants of the coefficient have insignificant intensities with respect to a certain 
threshold. Thus, due to this relationship, a chain of insignificant coefficients may be 

15 indicated by a single code, a technique that compresses the size of the data that indicates 
the original image. As an example, if the wavelet coefficient for the location 24 is a 
zerotree root, then the wavelet coefficients for the locations 20, 22 and 24 are 
insignificant and may be denoted by a single code. 

The coding of each decomposition level typically includes two passes: a dominant 

20 pass to determine a dominant list of wavelet coefficients that have not been evaluated for 
significance and a subordinate pass to determine a subordinate list of wavelet coefficients 
that have been determined to be significant. During the subordinate pass, a threshold 
may be calculated for each subband and used to evaluate whether coefficients of the 
subband are insignificant or significant. Unfortunately, due to the computational 

25 complexity, the above-described compression technique may be too slow for some 
applications, such as an interactive video compression application, for example. 

Thus, there is a continuing need for an arrangement that addresses one or more of 
the above-stated problems. 
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SUMMARY 

In one embodiment, a method includes providing wavelet coefficients that 
indicate an image and representing each wavelet coefficient as a collection of ordered 
bits. The bits of each order are coded to indicate zerotree roots that are associated with 
5 the order. 

Advantages and other features of the invention will become apparent from the 
following description, drawing and claims. 



BRIEF DESCRIPTION OF THE DRAWING 
10 Fig. 1 is an illustration of the hierarchical order of subbands produced by wavelet 

transformations. 

Fig. 2 is a schematic diagram of a computer system according to an embodiment 
of the invention. 

Fig. 3 is an illustration of a scanning path to determine zerotree roots according to 
15 an embodiment of the invention. 

Fig. 4 is an illustration of the organization of a wavelet coefficient matrix 
according to an embodiment of the invention. 

Fig. 5 is an illustration of a scanning path for a wavelet coefficient matrix. 
Fig. 6 is an illustration of a path that is traversed to locate zerotree roots. 
20 Fig. 7 is a flow chart illustrating the execution of a program to encode wavelet 

coefficients according to an embodiment of the invention. 

DETAILED DESCRIPTION 
Referring to Fig. 2, an embodiment 1 19 of a compression program in accordance 
25 with the invention may cause a processor 1 12 to encode wavelet coefficients in a bit- wise 
fashion. In this manner, instead of classifying the wavelet coefficients (as zerotree roots 
or isolated zeros, as examples), the processor 1 12 may produce codes to classify the bits 
of the wavelet coefficients. For example, in some embodiments, the processor 1 12 may 
classify a particular bit as being either a zerotree root, an isolated zero, a positive node or 
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a negative node. Unlike conventional zerotree coding schemes, thresholds are not 
computed to identify insignificant values, as the "0" bit is treated as being insignificant 
and the "-1" and "1" bits are treated as being significant. 

In this manner, the processor 112 may generate one of the following codes to 
5 classify a particular bit: a "P" code to indicate a positive node if the bit indicates a "1"; 
an "N" code to indicate a negative node if the bit indicates a an "R" code to indicate 
that a "0" bit is a zerotree root ; and an "IZ" code to indicate that a "0" bit is an isolated 
zero. In some embodiments, a particular bit is classified as a negative node only if the 
bit is the most significant nonzero bit and the bit indicates a "-1." For example, for a 

10 coefficient of "-3" that is represented by the three bits "-Oil," the processor 1 12 
generates an N code to represent the middle bit. However, for this example, the 
processor 1 12 generates a P code to represent the least significant bit. 

For purposes of providing the wavelet coefficients, the processor 112 may, via 
wavelet transformations, decompose coefficients that represent pixel intensities of an 

15 original image. These wavelet coefficients, in turn, form subbands that are located in 
multiple decomposition levels. To classify the bits, the processor 1 12, in some 
embodiments, may execute the program 1 19 to process the bits based on their associated 
bit position, or order. In this manner, the bits of each bit order form a hierarchical tree 
that the processor 1 12 may traverse to classify each of the bits of the tree as being either 

20 a zerotree root, an isolated zero, a negative node or a positive node. Thus, as an example, 
the most significant bits of the wavelet coefficients(this bit may also be zero) are 
associated with one hierarchical tree (and one bit order), and the next most significant 
bits are associated with another hierarchical tree (and another bit order). 

For example, if the absolute maximum wavelet coefficient is represented by three 

25 bits (as an example), then all of the wavelet coefficients may be represented by three bits. 
Therefore, for this example, three hierarchical trees are formed. In this manner, the 
processor 1 12 produces a code for each bit based on its indicated value (i.e., "-1," "0," or 
"1") and possibly (if the bit indicates a "0") its position in the associated hierarchical 
tree. 



5 




In some embodiments, the processor 1 12 indicates the P, N, IZ and R codes via a 
bit stream that progressively indicates a more refined (i.e., a higher resolution) version of 
the original image over time. For example, the processor 112 may use the bits "00" to 
indicate the "P" code, the bits "01" to indicate the "N" code, the bits "10" to indicate the 
5 "R" code and the bits "11" to indicate the IZ code. Other coding schemes are possible. 
The progressive nature of the bit stream is attributable to the order in which the processor 
1 12 processes the bit orders. For example, in some embodiments, the processor 112 may 
process the bit orders in a most significant first fashion. Therefore, the processor 112 
may initially produce code all the bits that have the highest bit order, then produce code 

10 for all of the bits that have the next highest bit order, etc. As a result of this progressing 
coding, the resultant bit stream may initially indicate a coarser version of the original 
image. However, more refinements to the image are indicated by the bit stream over 
time, as the processor 1 12 produces the codes for the bits having the lower bit orders. 
Thus, in some embodiments, the resolution of the image that is indicated by the bit 

15 stream improves over time, a feature that may be desirable for bandwidth-limited 

systems. As a result, a decrease in resolution of the reconstructed image may be traded 
for a decrease in communication bandwidth. 

Referring to Fig. 3, in some embodiments, the processor 1 12 process the bits of 
each order in a predefined sequence. For example, for a particular bit order, the 

20 processor 112 may begin with the highest decomposition level and produce codes for the 
bits of the highest decomposition level before proceeding to produce codes for the bits of 
the next highest decomposition level. The processor 112 produces code(s) for the bit(s) 
of the LL subband and, then for each decomposition level, produces code(s) for the bit(s) 
of the LH subband, subsequently, produces code(s) for the bit(s) of the HL subband and 

25 lastly, produces code(s) the bit(s) of the HH subband. 

As an example, the wavelet coefficients produced by a two level decomposition 
may be arranged in a matrix 40 that is depicted in Fig. 4. In this manner, the matrix 40 
may be viewed as being subdivided into four quadrants 30a, 30b, 30c and 30d. The 
upper right 30b, lower left 30c and lower right 30d quadrants includes the coefficients for 
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the LH, HL and HH subband images, respectively, of the first decomposition level The 
coefficients for the LL, LH, HL and HH subband images of the second decomposition 
level are located in the upper left 32a, upper right 32b, lower left 32c and lower right 32d 
quadrants of the upper left quadrant 30a. The coefficients produced by further 
decomposition may be arranged in a similar manner. For example, for a third level of 
decomposition, the upper left quadrant 32a includes the wavelet coefficients of the LL, 
LH, HL and HH subbands of the third decomposition level. 

If the coefficient matrix that indicates the pixel intensities for the original image 
is a 4X4 matrix, then the matrix 40 may be of the form that is depicted in Fig. 5. In this 
manner, the LL, LH, HL and HH subband images of the second decomposition level each 
have one coefficient, represented by "A" (for the LL subband image), "B" (for the LH 
subband image), "C" (for the HL subband image) and "D" (for the HH subband image), 
respectively. As depicted in Fig. 5, for the first decomposition level, the coefficients for 
the LH, HL and HH subband images are represented by the following respective 
matrices: 



l F3 



E, E 2 
E, E 4 



F, F 2 
L F 3 F J 



G, G 2 
G, G 4 



It is noted that each coefficient of the second decomposition level (except A), is 
20 associated with at least four coefficients of the first decomposition level, i.e., each 

coefficient of the first decomposition level has at least four descendant coefficients in the 
second decomposition level. Therefore, each bit in the first decomposition level has at 
least four descendent coefficients in the second decomposition level. 

For each bit order, the processor 112 may process the bits in the scanning 
25 sequence described above. If a particular bit indicates a "1" or a "-1," then the processor 
1 12 generates the P or N code and proceeds to process the next bit in the scanning 
sequence. However, if a particular bit indicates a "0," then the processor 1 12 may trace 
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the bit through its descendants to determine if the bit is an isolated zero or a zerotree root. 
The coefficients in the LL subband are simply entropy encoded. 

As an example, to produce the code for the least significant bit (called D(l)) of 
the D coefficient (located in the HH subband of the second decomposition level), the 
processor 112 determines whether the D(l) bit indicates a "0." If so, the processor 1 12 
evaluates the descendant bits G^l), G 2 (l), G 3 (l) and G 4 (l) of the subband HH of the first 
decomposition level in search of a "1" or "-1," as indicated in Fig. 6. If one of these bits 
indicates a "1" or "-1," then the D(l) bit is an isolated zero. Otherwise the D(l) bit is a 
zerotree root. 

As a numeric example, a 4X4 coefficient matrix that indicates pixel intensities for 
an image may undergo a two level decomposition to form the following matrix: 

'4 112" 
-2 0 0 1 

0 3 0 0 

0 1 0 0_ 

Because the maximum absolute value is "4," three bits may be used to represent the 
coefficients, as depicted in the following matrix: 



100 


001 


001 


010 


-010 


000 


000 


001 


000 


011 


000 


000 


000 


001 


000 


000 



Therefore, the processor 1 12 begins the encoding by generating codes for the third order 
bits (i.e., the most significant bits, which may be zero also) of the coefficients. More 
particularly, to generate the codes for the third order bits, the processor 1 12 follows the 
path 28 (see Fig. 5) and produces the appropriate code for the third bit of each coefficient 
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along the path 28. If a particular bit indicates a "0," then the processor 112 evaluates the 
descendents of the bit to find isolated zeros and zeroroots. The coding of the third order 
bits by the processor 1 12 produces the following codes (listed in the order of production): 
P,R,R,R. Subsequently, the processor 1 12 produces the codes for the second order bits 
5 (listed in order of production): IZ,IZ,N,R,IZ,P,IZ,IZ,IZ,P,IZ,IZ. Lastly, the processor 
1 12 produces the codes for the first order bits (listed in order of production): 
IZ,P 5 IZ,R 5 P 5 IZJZ,P 5 IZ,P,IZ,P. As described above, the processor 1 12 may indicate the 
codes via a two bit coding scheme and transmit the codes as produced via a bit stream. 
As an example, another processor 200 (see Fig. 2) may use the bit stream to 
10 reconstruct the coefficient matrix that indicates the pixel in intensities of the original 
image in the following manner. Before the decoding begins, the processor 200 first 
receives an indication from the processor 112 that three levels of coding (i.e., one level 
t3 for each bit order) have been used. After obtaining this information, the processor 200 

Uj may reconstruct the original coefficient matrix using the codes in the order that the codes 

15 are produced. More particularly, the processor 200 may use the codes produced by the 
;*f coding of the bits of the third bit order (i.e., the first level of coding) to produce the 

m following matrix: 



"100 


000 


000 


000 


000 


000 


000 


000 


000 


000 


000 


000 


000 


000 


000 


000 



20 

The processor 200 may use this matrix to reconstruct a coarse version (i.e., a lower 
resolution version) of the original image. However, if a more refined version is desired, 
the processor 200 may use the codes that are produced by the coding of the second bit 
order (i.e., the second level of coding) to produce the following matrix: 

25 
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100 000 000 000 

-010 000 000 000 

000 010 000 000 

000 000 000 000_ 

Finally, if the processor 200 uses the codes that are produced by the coding of the bits of 
the first order (i.e., the third level of coding), the processor 200 produces the original 
matrix of decomposed wavelet coefficients. 

Referring to Fig. 7, to summarize, the compression program 1 19, when executed 
by the processor 1 12 may cause the processor 1 12 to perform the following procedure to 
produce the above-described coding. First, the processor 1 12 may express (block 72) a 
matrix of decomposed coefficients in a signed binary representation. Next, the processor 
1 12 may determine (block 74) the number of digits that are needed to represent the 
absolute value of the maximum wavelet coefficient. This processor 1 12 uses a variable 
(called n) that indicates the current bit order being processed by the processor 112. In 
this manner, the processor 112 uses a software loop to process the bits, one bit order at a 
time. To accomplish this, the processor 112 produces codes (block 76) for the bits of the 
current bit order the using the techniques described above. Subsequently, the processor 
1 12 determines (diamond 78) whether the rate of transmitted bits may exceed a 
predetermined bit rate. If so, the processor 1 12 terminates the coding for the current 
image to comply with the predetermined bit rate. Otherwise, the processor 112 
determines (diamond 80) if all bit orders have been processed, i.e., the processor 1 12 
determines if n equals "1." If not, the processor 1 12 decrements (block 75) the order that 
is indicated by the n variable by one and proceeds to block 76 to traverse the loop another 
time to produce codes for the bits of another bit order. Otherwise, the coding is 
complete. 

Referring back to Fig. 2, in some embodiments, the processor 1 12 may be part of 
a computer system 100. The computer system 100 may include a bridge, or memory hub 
1 16, and the processor 1 12 and the memory hub 1 16 may be coupled to a host bus 1 14. 
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The memory hub 116 may provide interfaces to couple the host bus 1 14, a memory bus 
129 and an Accelerated Graphics Port (AGP) bus 1 1 1 together. The AGP is described in 
detail in the Accelerated Graphics Port Interface Specification, Revision 1.0, published 
on July 31, 1996, by Intel Corporation of Santa Clara, California. A system memory 118 
may be coupled to the memory bus 129 and store the compression program 119. As 
described above, the compression program 119, when executed by the processor 1 12, 
may cause the processor 1 12 to provide wavelet coefficients that indicate an image and 
represent each wavelet coefficient as a collection of ordered bits. The processor 1 12 
codes the bits of each order to indicate zerotree roots that are associated with the order. 

Among other features of the computer system 100, a display controller 113 (that 
controls the display 1 14) may be coupled to the AGP bus 1 1. A hub communication link 
1 15 may couple the memory hub 1 16 to another bridge circuit, or input/output (I/O) hub 
120. In some embodiments, the I/O hub 120 includes interfaces to an I/O expansion bus 
125 and a Peripheral Component interconnect (PCI) bus 121. The PCI Specification is 
available from The PCI Special Interest Group, Portland, Oregon 97214. 

A modem 140 may be coupled to the PCI bus 121 to a telephone line 142. In this 
manner, the modem 140 may provide an interface that permits the bit stream that is 
produced by the processor 1 12 to be communicated to the processor 200. The I/O hub 
120 may also include interfaces to a hard disk drive 132 and a CD-ROM drive 133, as 
examples. An I/O controller 117 may be coupled to the I/O expansion bus 125 and 
receive input data from a keyboard 124 and a mouse 126, as examples. The I/O 
controller 117 may also control operations of a floppy disk drive 122. Copies of the 
program 119 may be stored on, as examples, the hard disk drive 132, a diskette or a CD- 
ROM, as just a few examples. 

In the context of this application, the phrase "computer system" may generally 
refer to a processor-based system and may include (but is not limited to) a graphics 
system, a desktop computer or a mobile computer (a laptop computer, for example), as 
just a few examples. The term "processor" may refer to, as examples, at least one 
microcontroller, X86 microprocessor, Advanced RISC Machine (ARM) microprocessor, 
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or Pentium-based microprocessor. The examples given above are not intended to be 
limiting, but rather, other types of computer systems and other types of processors may 
be included in embodiments of the invention. 

Other embodiments are within the scope of the following claims. For example, 
the matrices of decomposed coefficients described above have one coefficient in each 
subband of the highest decomposition level. However, this arrangement is for purposes 
of simplifying the discussion of the coding. Therefore, each subband of the highest 
decomposition level may have multiple coefficients, and the above-described techniques 
may be applied to code the bits associated with these coefficients. In some embodiments, 
the processor 1 12 may code all of the bits of each order in parallel. In this manner, the 
coding of the bits of each bit order may be performed by the processor's execution of a 
separate thread. Other arrangements are possible. 

While the invention has been disclosed with respect to a limited number of 
embodiments, those skilled in the art, having the benefit of this disclosure, will 
appreciate numerous modifications and variations therefrom. It is intended that the 
appended claims cover all such modifications and variations as fall within the true spirit 
and scope of the invention. 
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